Guidelines for accurate genotyping of SARS-CoV-2 using amplicon-based sequencing of clinical samples (preprint)

BackgroundSARS-CoV-2 genotyping has been instrumental to monitor virus evolution and transmission during the pandemic. The reliability of the information extracted from the genotyping efforts depends on a number of aspects, including the quality of the input material, applied technology and potential laboratory-specific biases. These variables must be monitored to ensure genotype reliability. The current lack of guidelines for SARS-CoV-2 genotyping leads to inclusion of error-containing genome sequences in studies of viral spread and evolution. ResultsWe used clinical samples and synthetic viral genomes to evaluate the impact of experimental factors, including viral load and sequencing depth, on correct sequence determination using an amplicon-based approach. We found that at least 1000 viral genomes are necessary to confidently detect variants in the genome at frequencies of 10% or higher. The broad applicability of our recommendations was validated in >200 clinical samples from six independent laboratories. The genotypes of clinical isolates with viral load above the recommended threshold cluster by sampling location and period. Our analysis also supports the rise in frequency of 20A.EU1 and 20A.EU2, two recently reported European strains whose dissemination was favoured by travelling during the summer 2020. ConclusionsWe present much-needed recommendations for reliable determination of SARS-CoV-2 genome sequence and demonstrate their broad applicability in a large cohort of clinical samples.

Prospective mapping of viral mutations that escape antibodies used to treat COVID-19 (preprint)

Antibodies are becoming a frontline therapy for SARS-CoV-2, but the risk of viral evolutionary escape remains unclear. Here we map how all mutations to SARS-CoV-2's receptor-binding domain (RBD) affect binding by the antibodies in Regeneron's REGN-COV2 cocktail and Eli Lilly's LY-CoV016. These complete maps uncover a single amino-acid mutation that fully escapes the REGN-COV2 cocktail, which consists of two antibodies targeting distinct structural epitopes. The maps also identify viral mutations that are selected in a persistently infected patient treated with REGN-COV2, as well as in lab viral escape selections. Finally, the maps reveal that mutations escaping each individual antibody are already present in circulating SARS-CoV-2 strains. Overall, these complete escape maps enable immediate interpretation of the consequences of mutations observed during viral surveillance.

Role of Long-range Allosteric Communication in Determining the Stability and Disassembly of SARS-COV-2 in Complex with ACE2 (preprint)

Severe acute respiratory syndrome (SARS) and novel coronavirus disease (COVID-19) are caused by two closely related beta-coronaviruses, SARS-CoV and SARS-CoV-2, respectively. The envelopes surrounding these viruses are decorated with spike proteins, whose receptor binding domains (RBDs) initiate invasion by binding to the human angiotensin-converting enzyme 2 (ACE2). Subtle changes at the interface with ACE2 seem to be responsible for the enhanced affinity for the receptor of the SARS-CoV-2 RBD compared to SARS-CoV RBD. Here, we use Elastic Network Models (ENMs) to study the response of the viral RBDs and ACE2 upon dissassembly of the complexes. We identify a dominant detachment mode, in which the RBD rotates away from the surface of ACE2, while the receptor undergoes a conformational transition which stretches the active-site cleft. Using the Structural Perturbation Method, we determine the network of residues, referred to as the Allostery Wiring Diagram (AWD), which drives the large-scale motion activated by the detachment of the complex. The AWD for SARS-CoV and SARS-CoV-2 are remarkably similar, showing a network that spans the interface of the complex and reaches the active site of ACE2, thus establishing an allosteric connection between RBD binding and receptor catalytic function. Informed in part by the AWD, we used Molecular Dynamics simulations to probe the effect of interfacial mutations in which SARS-CoV-2 residues are replaced by their SARS-CoV counterparts. We focused on a conserved glycine (G502 in SARS-CoV-2, G488 in SARS-CoV) because it belongs to a region that initiates the dissociation of the complex along the dominant detachment mode, and is prominent in the AWD. Molecular Dynamics simulations of SARS-CoV-2 wild-type and G502P mutant show that the affinity for the human receptor of the mutant is drastically diminished. Our results suggest that in addition to residues that are in direct contact with the interface those involved in long range allosteric communication are also a determinant of the stability of the RBD-ACE2 complex.